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Exponential  Fourier  Densities  on  S /± 1 
and 

Optimal  Estimation  for  Axial  Processes 
by 

James  Ting-Ho  Lo  and  Linda  R.  Eshleman 
Mathematics  Research  Report  No.  77-1 
January,  1977 

ABSTRACT 

Several  models,  which  are  believed  to  be  generic  for  the  estimation  of 
discrete-time  axial  processes,  are  proposed.  By  introducing  axial  exponential 
Fourier  densities  and  axial  exponential  trigonometric  densities,  finite-dimen- 
sional recursive  schemes  are  obtained  for  updating  the  conditional  density 
functions.  The  underlying  idea  is  the  closure  properties  under  the  Bayes  rule 
of  the  various  combinations  of  these  exponential  densities. 

An  estimation  error  criterion,  which  is  compatible  with  a Riemannian  metric, 
is  introduced.  It  is  shown  that  the  corresponding  optimal  axial  estimates  can 
be  easily  computed  from  the  conditional  probability  distributions. 
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I.  INTRODUCTION 


In  this  paper  we  consider  the  problem  of  estimating  axes  in  three- 
dimensional  space.  An  axis  or  axial  vector  is  distinguished  from  a polar  vector 
in  that  the  former  is  invariant  under  inversion.  Such  axes  occur  in  many 
diverse  areas  including  the  following:  geophysical  fluid  dynamics  to  estimate 
the  vorticity  of  a flow,  paleoraagnetism  to  estimate  a magnetic  field,  crystal- 
lography to  estimate  the  optic  axis  of  a crystal,  geology  to  estimate  the 
direction  of  a normal  to  the  axis  of  a fold  in  a layer  of  rock,  and  quantum 
mechanics  to  estimate  the  axis  of  rotation  of  a rigid  body  rotation. 

Using  densities  of  the  form  exp  f where  f is  a linear  combination  of 
axially  symmetric  spherical  harmonics,  estimation  problems  which  arise  by 
examining  various  possible  ways  of  obtaining  a displacement  of  an  axis  will  be 

solved  in  this  paper.  Although  the  state  space  under  consideration  is  home- 

2 2 
omorphic  to  a hemisphere  of  S , the  results  in  [1]  for  estimation  on  S can- 
not be  applied  for  several  important  reasons:  the  displacements  defined  in 

that  paper  may  result  in  a given  point  being  displaced  to  a non-antipodal 
point  in  the  opposite  hemisphere,  the  densities  discussed  in  [1]  were  not,  in 

2 

general,  axially  symmetric,  and  the  error  criterion  used  for  S is  undesirable 
since  it  would  result  in  a rejection  of  the  antipode  of  the  optimal  estimate. 

Using  the  various  displacements  and  conditional  densities  obtained  in 
this  paper,  detection  for  axial  processes  would  be  described  by  procedures 
similar  to  those  used  for  S and  S0(3)  which  have  been  discussed  in  [1]  and  [2J. 

II.  PRELIMINARIES  AND  AXIAL  EXPONENTIAL  FOURIER  DENSITIES 

The  set  of  all  axes  centered  at  a fixed  point  can  be  identified 
in  a one-to-one  manner  with  the  points  on  a hemisphere  llS(p)  of  the  unit 
sphere  S with  pole  at  p.  Hence  the  state  space  for  the  axial  processes  to  be 
discussed  can  be  considered  to  be  the  factor  space  obtained  by  partitioning 


1 


2 2 
S into  equivalence  classes  such  that  two  different  points  on  S belong  to  the 

same  equivalence  class  if  and  only  if  they  are  antipodal  points.  This  factor 

2 

space,  which  is  projective  two-space,  will  be  designated  by  the  symbol  S /± 1 
throughout  this  paper. 

Since  the  spherical  coordinates  (0,4>),  (ir-9 , 4>+tt)  , and  (tt+0,4>)  will  all 

represent  tne  same  axis,  any  probability  density  p(0,<t>)  on  S2/±l  should  have 

2 

the  property  that  it  is  antipodally  symmetric  on  S , i.e. 

p(9»<t>)  = p(n-6,<t>+iT)  = p (tt+0 , 4>)  . (1) 

An  example  of  such  an  antipodally  symmetric  density  that  has  been  previously 
studied  is  the  Bingham  distribution  [3]  which  is  a generalization  of  the 
Dimroth-Watson  distribution  of  [4]  and  [5].  It  can  be  expressed  as  a density 
containing  spherical  harmonics  of  the  form 

2 

c exp  j a Y.  (9,<<>)  , 

m=-2  m 2’m 

where  c is  a normalizing  constant. 

2 

Clearly,  given  any  density  p^(0,<f>)  defined  on  S , a corresponding  density 

2 

p ( 0 , 4>)  can  be  obtained  on  S /±1  regarded  as  the  upper  hemisphere 
HS(n)  = [0,tt/2]  x [0,2tt] 

by  defining 

p(0,<|>)  ” P^Q.'f1)  + p^TT-O.tH77)  » (9»4>)  f HS(n) 

However,  having  seen  in  [1]  that  exponential  Fourier  densities  were  most  use- 
ful in  determining  conditional  densities  that  are  needed  for  estimation,  we 
shall  instead  consider  those  densities  of  the  form 

2 


M)>. 


N £ 

p(6,4>)  = exp  l l a Y (0,<j>) 
£=0  m=-£  m m 


for  which  (1)  holds.  For  9 in  the  interval  [0,it]  the  spherical  harmonics 


Y.  (0,(1))  are  defined  by 
£m 


W0^  " ^ 


f ( £-m) ! 1 ~ 
|_(£+m)  ! j 


P™(cos  0)eim<^ 


where  P (x)  is  an  associated  Legendre  function. 

Xj 

We  extend  this  definition  for  0 outside  this  interval  by  defining  P™(cos  0) 
in  such  a way  that  the  values  of  the  spherical  harmonics  remain  independent  of 
the  choice  of  spherical  coordinates.  In  particular,  for  -tt  < 0 < 0,  (2)  becomes 


V-M+’> ' [&st]1/2  p"(c°s  <-9))e 


so  that  in  order  for  the  identity 


■ Ve’*> 


to  hold  we  must  define 


P™(cos  0)  = (-1)P1  P™(cos  (-0))  for  -ir  < 0 < 0 
Similarly,  we  define 

p“(cos  0)  = (-l)m  P™(cos  ( 2tt— 0 ) ) for  tt  < 0 < 2n 

Xj  X, 

Suppose  we  now  consider  the  point  (0-Ht ,<J>)  which  is  the  antipodal  point 


of  (0,<|>)  where  0 < 0 < i.  We  have 
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Y._(0+M)  = (-i)m  I^glj172  p”(cos  (0-Hr))elra,> 


-1/2 


P^ (cos  (TT-e))e 


m / « s 1 m4> 


P^(-cos  0)e 


(-D 


nrt-2 


f (2-m)  !~] 
[_(fi.-Hn)  !j 


1/2 


p”(cos  0)eim^ 


■ (-1’  V9-*1  • 


The  fourth  equality  follows  from  the  property  of  associated  Legendre  polynomials 
that 


P™(-cos  0)  = (-l)®**  P*(cos  0)  for  0 £ [0,TT] 


Since  the  spherical  harmonics  are  independent  of  the  choice  of  coordinates 
we  can  conclude  that  any  spherical  harmonic  assumes  the  same  value  at  ant  i- 


podal  points  if  and  only  if  it  is  of  even  order  and  hence  we  have  the  following 

,2 


lenima  which  characterizes  all  KFD's  for  which  antipodal  points  of  S have 
the  same  density: 


Lemma  1:  Let  p(0,<{>)  = exp  f(0,<J>)  where 


- I ! ve.*> 

2=0  m=-2 


Relation  (1)  holds  if  and  only  if  a„  =0  whenever  2 is  an  odd  integer. 


Thus  we  define  an  axial  exponential  Fourier  density  of  order  2H  on  S /±1, 
denoted  by  AEFD(2N) , to  be  a density  of  the  form 


N 22 


p(0,4>)  - exp  l l a 


£-0  m— 22, 


22, m W6’^ 


(3) 


— J 
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In  the  special  case  where  the  density  is  isotropic  so  that  p(0,<t>)  is  inde- 
pendent of  <(>  (3)  is  replaced  by 

N N 

P(6)  = exp  l a2Jl  Y2j,f0^e’^)  = exp  £ a2l  P2SL(c°S  0) 

It  should  be  noted  that  any  twice  continuously  differentiable  density  can 

2 

be  approximated  by  an  AEFD  in  the  space  of  square  integrable  functions  on  S /±1 
as  a consequence  of  Theorem  1 in  [1]. 

In  the  determination  of  conditional  densities  for  the  several  estimation 
models  to  be  discussed,  it  will  frequently  be  necessary  to  transform  densities 
of  the  form  (3)  into  expressions  containing  the  trigonometric  functions 
exp[i(m0+n<}>)  ] . By  Lemma  3 of  [1]  a spherical  harmonic  of  even  order  can  be 
expressed  as  a linear  combination  of  trigonometric  functions  of  the  form 

N 2N 

exp  l l b exp  [i(2j0+k<f>)  ] (4) 

J—  N k=-2N 

which  we  will  refer  to  as  an  axial  exponential  trigonometric  function  of  order 
2N,  denoted  by  AETF(2N)  to  differentiate  it  from  an  ETF(2N)  of  [1]  where  the 
coefficient  of  0 can  be  an  odd  integer.  It  should  be  observed  that  the  use  of 
the  word  "axial"  is  a misnomer  since  (1)  does  not  hold  for  densities  of  the 
form  (4) . 

Our  choice  of  an  error  criterion  for  optimal  estimation  of  a random  axis 

is  motivated  by  the  desire  that  antipodal  points  be  indistinguishable.  If  x 

2 

and  y are  the  column  vectors  of  direction  cosines  of  two  points  on  S and 
p(x,y)  denotes  the  principal  angle  between  them,  then  we  shall  use  the  measure 


1 - cos  2p(x,y)  ■ 2 [1  - cos2  p(x,y)]  = 2[1  - (x'y)2l 


NBRR! 


fore,  if  x represents  a random  direction  cosine  vector,  the  axial  estimation 
problem  will  be  to  find  its  estimate  y which  minimizes  the  criterion 


J(x,y)  = 2 E [1  - (x'y)2]  = 2[1  - y'E(xx’)y]  (5) 

, 2 

where  the  estimation  is  over  an  arbitrary  hemisphere  of  S . 

Since  E(xx’)  is  a positive  definite  3x3  symmetric  matrix,  it  is  known 
that,  just  as  in  the  S0(3)  case,  the  quadratic  form  y'E(xx’)y  is  maximized 
when  y is  the  unit  eigenvector  associated  with  the  largest  eigenvalue  of  E(xx'). 

Hence 

min  J(x,y)  = 1 - x'E(xx')x  = 1 - A 

y 

' 

where  A is  the  maximum  eigenvalue  of  the  matrix  E(xx')  and  2 is  the  unit 


eigenvector  associated  with  A. 

Once  x'  = [2^  2^  is  calculated,  the  spherical  coordinates  (0,<J>) 

of  the  optimal  estimate  can  be  determined  from  the  relations 


2 

III.  DISPLACEMENTS  ON  S / ±1 

2 2 

In  order  to  use  a hemisphere  of  S as  a representation  for  S /±1  in  the 

estimation  of  axial  processes,  it  is  necessary  to  use  displacements  o which 

2 — 

have  the  property  that  for  any  xf  S /±1,  the  displaced  points  x O y and  x o y 

2 — 

are  either  identical  or  antipodal  points  on  S , x being  the  antipode  of  x. 


While  the  displacements  used  for  estimation  on  S in  [1]  do  not  have  this 

property,  it  will  be  seen  that  variations  of  these  displacements  which  have 

this  property  can  be  defined  which  result  in  usable  estimation  models. 

2 

Let  us  consider  the  first  displacement  on  the  sphere  S with  center  0 
given  in  [1]  of  a point  s (6 » 4>)  by  an  ordered  pair  v(a,B).  Assuming  that  (0,<f>) 
and  (a, 3)  are  both  in  [0,it)  x [0,2tt),  let  n denote  the  point  (0,0,1)  and  let 
ns  be  the  smaller  great  circle  arc  between  n and  s.  There  is  a unique  semi- 

/ — v 

circle  T with  center  at  0 having  s as  an  endpoint  such  that  the  angle  from  ns 

to  T has  measure  3,  using  the  right-hand  screw  rule  in  the  direction  Os.  The 
point  m which  is  the  endpoint  of  the  arc  sm  on  T of  length  a is  shown  in 

Figure  1 and  its  coordinates  (A,y)  are  related  to  those  of  s and  v by  the  equations 


cos  a = cos  A cos  0 + sin  A sin  0 cos  (y-4>) 

sin  a cos  3 = cos  A sin  0 + sin  A cos  0 sin 

sin  a sin  3 = sin  A sin 


By  considering  the  special  case  where  3 = 0 it  is  seen  that  the  spherical 

coordinates  of  m are  (0-a,<{>);  but  the  displacement  of  the  signal's  antipode 

_ - 2 
s(tt-0,4>+tt)  is  (TT-0-a,4>+^)  which  does  not  coincide  with  either  m or  m on  S . 


Figure  1 


This  difficulty  can  be  resolved  by  making  the  restriction  that  a fixed  hemi- 
sphere be  used  as  the  state  space.  For  convenience  we  choose  the  upper  hemi- 
sphere and  require  that  its  points  have  spherical  coordinates  (6,<p)  such  that 
0 < 0 < it/2  and  0 < <f>  < 2tt.  Throughout  the  remainder  of  this  paper  we  will 
identify  S^/±l  with  this  set 


S /±1  = [0,tt/2]  x [0, 2tt)  . 


If  s(0,<(>)  and  v(a,B)  belong  to  S^/il  let  n(A,y)  be  the  displaced  point 
where  (A,y)  is  in  [ — tt , tt)  x [-2h,2tt).  Now  we  define  m(A,y)  as  being  a point  of 


S /±1  with  coordinates 


<A,y)  = 


(n+A.y  mod  2tt) 
(-A,(u+ti)  mod  2tt) 
(A,y  mod  2tt) 
(tT-A,(y+TT)  mod  2rr) 


if  A £ -it/ 2 
if  — tt/ 2 < A £ 0 
if  0 < A £ tt/2 
if  A > 71  / 2 


and  we  indicate  this  relation  of  m to  s and  v by 


■ £ 


m(A,y)  = s(0,4>)  ^a^  v(ot,B) 

Geometrically,  it  is  seen  that  ^a^  is  identical  to  the  spherical  displace- 
ment whenever  m is  in  HS(n),  while  m is  the  antipode  of  m otherwise.  The  require- 
ment that  the  first  spherical  coordinate  be  between  0 and  tt/2  is  made  to  avoid 
some  unnecessary  difficulties  in  the  calculation  of  conditional  densities  to 
be  discussed  later. 

A second  displacement  © arises  by  increasing  each  coordinate  of  the 

2 

position  of  the  signal.  Let  Sq^q^q)  be  in  S ^±1  and  suPPose  ls  a 

pair  of  angles  satisfying  the  condition 

(Z,0  ( [-TT/2.TT/2)  X [0, 2tt)  . 


8 


Now  set  6^  = 0g  + £ and  <p^  = 4>q  + Let  be  the  point  of  S / ±1 


with  coordinates 


r(-e1,($1 

(0l»4>i)  = “I  (Qj,^  m 


-Hr)  mod  2tt) 
mod  2tt) 


(tt-0^,  (c^+tt)  mod  2tt) 


if  -tt/2  < 0^  < 0 
if  0 < ^ < tt/2 
if  tt/2  < 0^  <_  Ti 


This  relationship  of  to  s^  is  denoted  by 


Si(0i,<!)i)  = so(6o,<t>o)  ® 


If  0 < £ < tt/2  then,  since  0 £ 0^  + E,  £ tt,  the  two  hypotheses 


H1 : (0i» = (0q+5>  (^q+O  mod  27T) 

H2:  (ex » = (tt-0o-C,  (<()0+^+tt)  mod  2tt) 


must  be  considered  in  the  determination  of  s^  if  s^  and  w are  known.  If, 
instead,  -tt/2  < E,  < 0,  these  hypotheses  are  replaced  by 


H1 : (01 , 4>1)  = (0q+£,  (4>q+C)  mod  2tt) 

H2:  ©v^V^10  mod  2tt^ 


since  -tt/2  < 0^  + £ _<  tt/2. 

It  should  be  observed  th  it  ^b^  identifies  s^  with  the  antipode  of  s^  when 
s^  is  not  in  HS(n)  and  is  similar  to  the  second  displacement  © on  S2  in  [1] 
when  s^  is  in  HS(n). 

The  final  displacement  © to  be  used  is  obtained  from  a rotation  of  tnree- 

2 

dimensional  space.  If  sQ  is  in  S / ±1  and  R is  a rotation  that  rotates  sQ  to 

-2  ~ 

the  position  s^  on  S , we  define  s^  to  be  s^  if  is  in  HS(n)  and  to  be  the 
antipode  of  s^  if  s^  is  not  in  HS(n)  and  write 

S1  “ s0  ® R • 


9 


It  should  be  noted  that  since  R will  rotate  the  antipode  s^  into  s^,  unlike 
the  previous  displacements,  it  is  actually  unnecessary  to  require  that  the 
fixed  hemisphere  (7)  be  used  as  the  state  space. 

2 

IV.  ESTIMATION  OF  AN  AXIS  ON  S /+1 

2 

We  now  consider  a signal  s,  represented  by  the  point  s(0,0)  on  S /±1, 

2 

which  we  want  to  estimate,  given  an  observation  m(X,y)  in  S /±1  and  using  the 
error  criterion  (5).  The  optimal  estimate  requires  the  calculation  of  the 
conditional  density 

p(s)m)  = c p(m|s)  p(s)  (8) 

according  to  Bayes'  rule.  Because  of  the  nature  of  the  definitions  of  the 

various  displacements  in  the  previous  section,  m may  be  the  antipode  of  the 

2 

measurement  of  some  displacement  on  S , and  consequently  care  must  be  taken  in 
determining  the  density  p(m|s). 


Estimation  Model  E(a) 

2 

Let  m be  a measurement  on  S /±1  of  an  axial  signal  s which  can  be  described 
by  the  relation 

i(X,y)  = s(0,0)  © v(a,B) 


mi. 


where  s and  v have  AEFD(2N)'s  exp  f (0,0)  and  exp  f (ot, 8),  respectively,  where 

s v 


N U 

f (9,0)  = l l a 

£=0  m=-2£ 


21, m Y2«,,m(0,<1)) 


and 


N 21 

f(a.B)  - l l b Y (a,B)  . 

V £=  0 m=-2£  2*”m 

2 

Consider  the  possible  positions  of  m(X,y)  on  S /± 1.  First,  suppose  m(X,y) 


is  in  the  hemisphere  HS(s)  whose  pole  is  at  the  point  s(0,0).  Since 
0 < a < tt/ 2 , it  follows  that  we  have 


p(m|s)  = exp  f^Ca.g)  = gv(A,p) 


(9) 


where  g (A,y)  is  obtained  by  using  the  spherical  trigonometric  relations  in  (6). 

On  the  other  hand,  if  m(A,b)  is  not  in  HS(s),  then  it  must  be  the  antipode 
of  a point  m(A,u)  not  in  HS(n)  where 

A = ir  - A , y = (m+tO  mod  2n 

But  this  means  that  there  is  a unique  directional  process  v(a,B)  such  that 
m(A,y)  = s(0,<f>)  ^a^  v(a,8) 

where  0 < a < tt/2  and  hence 

p(m | s)  = p (m | s)  = exp  f~(a,8)  = g^(A,y)  . (10) 

If  gv(A,y)  = g~(A,jj)  then  p(m|s)  could  be  determined  in  the  same  manner 
as  in  model  SF-1  of  [1].  The  following  lemma  gives  the  important  result  that 
any  AEFD  in  the  variables  (a, 8)  has  this  property: 

Lemma  2:  If  Z„  (A,y)  denotes  the  function  obtained  by  substituting  the 

3cm 

relations  (6)  into  the  spherical  harmonic  Y^a.B)  then 

wx-u)  ■ (-i)l1  V-x>,J+">  • 

Proof:  By  examining  the  relations  (6)  we  see  that  the  statement  holds  for  the 

spherical  harmonics  of  orders  zero  and  one: 

Z00(X‘M)  = 1 

Z,,(A,u)  - - — (sin  a cos  8 + i sin  a sin  8) 

11  n 

Z1Q(A,w)  = cos  a 

Z.  (A,u)  = — (sin  a cos  8 - i sin  a sin  8) 

1,-1  Jl 
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Now  assume  the  statement  holds  for  all  i <_  L.  From  the  recursive  relations 
for  spherical  harmonics  we  have 


ZL+1  m(X,u)  = [ (L+m+1) (L-m+1) ] 1/2[(2L+l)Z10(X,u)ZLpi(X,y) 


- /(L-m) (L+m)  ZT  (X,u)] 
L-l,m 


when  m = -L,-L+l, . . . , L-1,L  and 

(X,p) 


Z1;L(X,p)  Zll(X,m) 


L+l.L+l 

If  L is  even,  it  is  evident  from  the  induction  hypotheses  that 


while  for  L odd  we  have 


L+l,m 


Hence  the  statement  holds  for  all  integers  Z. 


□ 


This  lemma  assures  us  that  (9)  and  (10)  have  the  same  value  since 
exp  f (a, 3)  is  an  AEFD  and  hence  without  any  loss  of  generality  we  can  assume 
that  m(X,y)  is  in  HS(s)  with  a and  3 defined  by  (6). 

By  Lemma  1 of  [1],  exp  fv(a,3)  is  an  exponential  trigonometric  function 


of  the  form 


2N 

exp  7 b exp  [i(m9+n<}>)  ] 
m, n=-2N 


when  transformed  by  (6)  into  a function  of  (0,4>)  so  that 


p(s|m)  = c exp  [fg(0,<|))+fv(ot,3)  ] 


is  also  an  exponential  trigonometric  function. 


We  remark  that  if  fv(a,0)  has  coefficients  with  the  special  property  that 
b2£,m  = (-1)  b2£,-m 

then  by  Lemma  3 of  [1],  exp  fv(a,0)  when  transformed  into  a function  of  (0,<f>) 
is  an  AETF(2N)  and  hence  p(s|m)  is  an  AETF(2N) . 

For  the  isotropic  noise  case  where 


N 

f (a,  0)  = £ 

v £=0 


P2Jl(cos  a) 


N 

I 


2£ 

l 


£=0  m=-2£ 


b2£  Y2£,m(A,M) 


Y2£,m(6-») 


it  is  observed  that  p(s|m)  is  itself  an  AEFD(2N) . 

A generalization  of  the  noise  density  to  include  the  situation  where  the 
density  depends  not  only  on  the  spherical  coordinates  (a,0)  but  also  on  the 
position  of  the  signal  s(0,<f>)  can  be  used.  Suppose 


f (a,  0) 
v 


0 2£ 

l l 

£=  0 m=-2£ 


where  each  bon  (8,4>)  is  a smooth  function  of  0 and  cj)  on  S / ±1  so  that  there 

2£,m 

is  an  integer  N(£,m)  and  coefficients  b(2£,m,2k,j)  such  that 


N(£,m)  2k 

J l b(2£,m, 2k, j ) Y ,(0,<t>) 
k=0  j=-2k  ZK,J 

approximates  b2^  m(0,<}>)  sufficiently  closely.  Since  Y2^  m(cx,8)  an^  ^2£  m^’^^ 
can  each  be  expressed  as  linear  combinations  of  terms  of  the  form  exp  (i(2j0+k<fi)  ] , 
their  product  and  hence  fv(a,0)  are  also  such  linear  combinations.  Therefore, 
the  noise  density  is  an  exponential  trigonometric  density  of  order  2M  where 

M = N_  + max  N(£,m) 


L3 


If  this  generalized  density  has  the  additional  property  that 


= (-1)  b90  * 

zx.,m 


then 


+ b22,m(8-»>Y2t,„(“-8)  ’ 

N(«,,m)  2k 

I l Y2k>j(0,(t)){b(2Jl,n.,2k,j)[Y9p  m(a,e)  + (-l)mY2p  _m(a,f 


k=0  j=-2k 


' 2£.  ,m 


24, -m' 


By  Lemma  3 of  [1],  we  can  thus  conclude  that  fv(a,3)  is  of  the  form 


0 11  N(£,m)  2k 

f (a, 3)  ‘III  l C(2£,m,2k,j)exp[i(2k6+j(),)]Y  (0,<t>) 

v 4-0  m=0  k=0  j=-2k  ’ 3 


and  hence  that  exp  f (a,  3)  is  an  AETF(2M)  in  (0 , <J>) . 

v 

2 

V.  ESTIMATION  OF  A TIME-VARYING  AXIAL  PROCESS  ON  S /±1 
Using  displacements  © and  © to  define  a sequence  of  signals  in 
conjunction  with  measurements  obtained  by  using  the  displacement  © as  well 
as  tne  Important  case  of  additive  white  noise,  we  now  describe  several  time- 
varying  estimation  models  tor  axial  processes. 


A.  Estimation  Model  E(b,a) 

Suppose  signal  and  measurement  processes  are  described  by 


WWW  = 8k((W  © vw* 
VW  = sk(ek’V  © Vvek> 


where  -it/2  < L < tt/2,  0 <_  < 2tt,  and  the  densities  of  sQ  and  v^  are  AEFD(2N)'s. 


Let  us  assume  that  0 £ £q  < tt/2  so  that 


p(Sl)  - p(s1|H1)P(H1)  + p(Sl|H2)P(H2) 
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lT  & 

r • 


and 


Let 


and 


f2n  ,t r/2-C 

P<»1>  ■ \0  J0  ><*o>  dso 


/2-rr  rTT/2 


p(»2) 


0 jtt/2-S 


P(80)  ds0 


gj^s.m)  = 


g2(s,m)  = < 


0 


1 if  0 > C 

0 otherwise 

1 if  0 >_  tt/2  - i 
0 otherwise 


so  that  if  p(sQ)  = p(sQ |m  ) = exp  f1(e0»(J>0)  then 

|-2tt  |-tt/2-C0 

P(Sl|m°)  = g1(s1,m0)exp  J Jq  P<*0>  ds0 


r2v  rTT/2 

+ g2(s1,m0)exp  f1(^-01-Co,4'1-V7r)Jo  Jv/2_£  P(S0} 


ds. 


Now  if  f j 


is  the  exponent  of  an  AEFD  so 


that  it  can  be  written  in  the  form 


then 


where 


N 2N  _ 

f _ (0 , 4>)  = y y a exp  [ i ( 2m0+n<)> ) ] , 

i _ «T  _ 1,1,1 

m=-N  n=-2N 


r N 2N 


p(Sl|m  ) - RQ1  exp 


l l amnexP  [-i(2mCQ-n^Q>  ]exp[i(2m0^+n4)1)  ] 


m=-N  n=-2N 
r N 2N 


+ R02  exp 


l l {(-l)a  exp[i(2m^n-nr  )]  exp[i(2m0  +n<})  ) ] } 

m-N  n-2N  "m’n 


f 2ti  ,t t/2-£0 

" 8i<8i»mn> 


'm 


P(sn)  ds 


R02  " S2(S] 


2tt  c7r/2 

0 Jtt/2-C, 


P(s0)  dsQ 


Since  pOn-Js^)  is  a density  of  the  form 


' 2N  2N  , 

exp  l £ bm  “PliW^)! 

m=-2N  n=-2N 

the  conditional  density  pCs-Jm1)  will  be  the  sum  of  two  exponential  functions: 

1 l f 2N  11 

p(s  |n>)  = c {R  exp  £ a^  exp  [iO^+n^)  ] 

I |_m,n=-2N 


+ Rno  exp  l a exp[i(m0  +ni}>  )] 
02  __  mn  i J- 

m,  n=-2N 


a°  exp[-i(m^n+n^  )]  + b* 
mn  0 0 mn 


if  m is  even 


if  m is  odd  , 


(-l)n  a0^  n exp[i(mC0-nC0) ] + b*n  if  m is  even 


if  m is  odd 


i 

A 

I 

I 


If  instead  -tt/2  <_  < 0,  then 


f2TT  (-tt/2  ["'  f-^Q 

P(H.)  = | p(s  ) ds  P(H  ) = P(s  ) dsQ 

1 J0  u ’ JoJo 


2TT  r-Zr 


so  that  if  we  let 
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g3(s,m)  = 


g4(s,m)  = 


if  0 < tt/2  + £ 


otherwise 


if  0 < -5 


otherwise 


o r27T  r12 

P(sj_  lm  ) = 83(8i»®o^  exp  ' j J p^s0^  dS0 

0 "f’0 

(2tt  rTT/2 

+ g4(srm0)  exp  f1(-ei-C0,<()1-C0)  * j J P(®0)  ds0 

0 “^0 


so  that 


1 r 2N  ii 

p(s  |m  ) = cx  RQ1  exp  £ aran  exp  [1(00^^)  ] 

(m,n=-2N 


+ Rno  exp  l a exptKrae.-HK})  ) ] j , 
Uz  u „„  mn  i J-  i 

m,n=-2N  | 


where 


a°  exp [-i(m£n+n£n) ] + b1 
ran  U U mn 


a°  n exp[i(m£  -nr  )]  + b1 
—in  }n  U U mn 


if  m is  even 


if  m is  odd  , 


if  m is  even 


if  m is  odd  , 


where  RQ1  * g3<81»m0)P(Hl) ’ and  RQ2  = ' The  computation  of 

p(s2|m2)  = p(s2|m1,H1)P(H1|m1)  + p(s2 |m1,H2>P(H2 lm1) 


17 


will  not  be  included  here  although  it  proceeds  in  a similar  manner.  We  observe 

k k 

that  for  any  k,  p(sk|m  ) will  be  a linear  combination  of  2 terms  of  the  form 


r 2N 


exp 


l exp  [i(m0. +n<(>,  ) ] 


m,n=-2N 


mn 


k'“Tk' 


k-1. 


and  it  will  be  necessary  to  determine  twice  as  many  parameters  as  for  p(skl|m  ) 


B.  Estimation  Model  E(c,a)  Using  Isotropic  Noise 

Suppose  the  signal  process  {s^}  is  described  by 


s—  = s’  © Rk 


k+1  k 


where  {R^}  is  a sequence  of  rotations  in  three-dimensional  space.  From  the 
definition  of  we  have,  using  direction  cosine  vectors  for  sQ  and  s^  and  the 
direction  cosine  matrix  representation  for  Rq 


so  = RoX  S1 


or 


so  = Ro1  S1 


where  s.^  is  the  antipode  of  s^  having  spherical  coordinates  ( 0^+tt  , tj^ ) . By 
[6,  p.150] 


2£ 


V»‘VV  ■ i , D <R0>™, 


n=-2£ 


for  either  case,  so  that  if 


N 21  _ 

P(sq)  ” exP  I 1 a 2 5J,  m Y2£  n/^0’^0^ 

0 £=0  m=-2£  X,,m  Z ,m  U 


then 


P(sL) 


N 2£ 

e*P  l l 

£-0  m--2£ 


' 2£ 

l 

n=-2£ 


0 _2£ . , 

3 2 £ , m°  R0  mn 


Y2£,m(<W 
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■ -rl— 


then 


p(V  = exp  J0  b2£  P22.(cOS  ak) 

N 21 


*\l\>  ■ “»  l J,>  Vm(W  Y2t..(W 

X/—  u m ^ X/ 


Therefore, 


N 2SL 


p(sk|mk)  = exp  I J an,mY2£,m(Xk‘yk) 
x,  0 m X/ 


where 


21 


a2SL,m  " ^ a2£,n  D (Rk-l)mn  + b2£  ^Jl.m^k’^ 

n=-2£ 


C.  Estimation  Model  E(b,+)  with  Additive  White  Noise 
Let  {s^}  be  a sequence  of  signals  such  that 


sk+l(0k+l,l*!k+l)  = sk(W  © wk(?k’V 

and  suppose  {m^}  is  a sequence  of  vectors  such  that 


mk  “ + vk 


where  {v^}  is  a sequence  of  independent  Gaussian  vectors  with  mean  zero  and 

2 

covariance  matrix  = Elv^v^l  an^  h is  a vector-valued  function  on  S /+ 1.  The 

2 

completeness  of  the  spherical  harmonics  of  even  order  on  S /+  1 allows  us  to 
consider  h(s)  to  be  approximately  of  the  form 


N 21 

l l h 

£-0  m—2*. 


2£,m  Y2£,m(6,<,>) 
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MfCgwr 


Now 


p(vk)  = (2TT)-p/2(det  Rk)-1/2exp[-  | ^ v£  vd  ] 


where  Rk  has  matrix  elements  Rj^  and  vector  vk  has  entries  v We  therefore 
have  as  the  conditional  density 

-17  #9  r N 2£ 

P(\|sk)  " (2lT)_P  <det  V"  exp  V ^ ^ Cka*m)Y2£,m(6k*<,)k) 

L 36=0  m=-236 


21  2SL 


11  4A/  4A/ 

+ , L J„  ,i„.Cka’,t'’,">”')Y2t,»<ek’V  Y2f,m(W 

36  j 36— U m ^36  m — Z36  . 


(ii) 


with 


1 ? i j 

ck  " ~ 2 . \ \ “k 

^ f J 


Cka’m)  = 1 . tRkj  [I\  h2£,m  + "4  h2£,m] 

1 » J i 


ck(£,  £‘,m,m*) 


i,j=l 


viJ  h1  i,3 

Xk  n2£,m  2i',m 


(11)  is  an  AETF(4N)  since  each  spherical  harmonic  of  even  order  is  an  AETF  and 
hence  so  is  the  product  of  two  such  spherical  harmonics. 

If  sQ  has  an  AEFD(2N) , then  using  (11)  the  conditional  density,  obtained 
in  the  same  manner  as  in  model  E(b,a),  will  be  a finite  linear  combination 
of  exponential  trigonometric  functions  of  order  4N. 


D.  Estimation  Model  E(c,+) 

As  our  final  model  consider  the  signal  and  measurement  processes  given  by 

WWW  = WV  © *k 

"k  " h(sk>  + vk 
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where  h and  v^  are  defined  as  for  model  E(b,+), 


If  p (m^ | s^)  given  by  (11)  can  be  written  as  an  AEFD(4U)  of  the  form 


2N  28, 


P^IV  = exp  l ^ b28,m  Y2£,m(ek’\) 

8=0  m=-28 


and  if 


2N  28  Q 

p(s0)  = exp  I l a2£>m  Yn>m(90,V  , 

x,=U  m=—  Z x. 


2N  2fL 

p(sk|»k)  - exp  I l a2a,,n  Y2H,m(W 
x,=U  m — zx. 


where 


'21,.  ' Jn  *£„  + b2*  WV"k> 


The  apparent  difficulty  in  requiring  that  (12)  holds  is  resolved  by  the 
following  lemma: 

Lemma  3:  The  product  of  two  spherical  harmonics  of  even  orders  is  a sum 

of  spherical  harmonics  of  even  orders. 

Proof:  The  product  of  two  spherical  harmonics  of  orders  28  and  28'  where 

8'  > 8 is  specified  by  [6,p.l65]: 


28+28 

Y2H,n/0,^)  Y2t',.<6’'H  ‘ LJLn 


where 


af  _ = (-1)  (2L+1) 

L,m+m 


28  28*  L \ /28  28*  L 


m -m-tn1  / \ 0 0 0 
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Now  the  proof  follows  as  an  immediate  consequence  of  Lemma  1. 

VI.  CONCLUSIONS 

Several  models  for  the  estimation  of  discrete-time  axial  processes  have 
been  proposed  and  have  been  found  to  be  solvable  by  introducing  axial  expon- 
ential Fourier  densities  which  employ  spherical  harmonics  of  even  orders. 

When  these  densities  are  expressed  as  trigonometric  functions  a closure  property 
for  the  conditional  densities  is  seen  to  exist  that  permits  their  computation 

by  updating  a finite  number  of  parameters. 

2 

The  three  displacements  defined  on  S /± 1 which  were  used  to  obtain  the 

estimation  models,  together  with  the  properties  of  spherical  harmonics  of 

even  order  yielded  conditional  densities  which  resembled  the  conditional 

2 

densities  arising  in  estimation  models  for  directional  processes  on  S that 

2 2 

appear  in  [1].  It  should  be  kept  in  mind,  however,  that  S and  S /±1  have 
inherent  geometrical  differences  as  is  evidenced  by  the  need  to  use  different 
error  criteria  on  the  spaces  involved. 

The  displacements  defined  for  S^/±l,  with  the  only  exception  being  , 

2 2 
required  that  S /±1  be  represented  by  a fixed  hemisphere  of  S . The  question 

arises  as  to  whether  or  not  it  is  possible  to  define  other  displacements  for 

which  this  requirement  can  be  removed. 
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